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Abstract 

Possibility of description of the glass transition in terms of critical dynamics considering 
a hierarchy of the intcrmodal relaxation time is shown. The generalized Vogel-Fulcher law 
for the system relaxation time is derived in terms of this approach. It is shown that the 
system satisfies the fluctuating-dissipative theorem in case of the absence of the intermodal 
relaxation time hierarchy. 

1 INTRODUCTION 

Despite of significant success in description of the dynamics of glass-transition, which they 
managed to attain by now, in particular in the framework of the mode coupling theory, the 
idea of description of glass-transition within the scaling theory of phase transitions stays as 
very attractive one. At first sight from the point of view of technical approaches such task 
does not present any complexity [1]. However, attempts of description of the phase transitions 
in disordered systems within the framework of the scaling approach faced serious difficulties 
connected with the necessity of a deviation from the framework of usual representations of the 
fluctuation theory of phase transitions [2] . 

It is possible that part of these difficulties is caused by using of the static fluctuation theory 
of phase transitions and by the complexity of averaging over configurations of the exchange 
integrals. In this work an attempt to consider this problem within dynamic approach, which 
is an alternative to the replica method and in which such difficulties are absent, is made. It 
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is assumed, what the ultrametric valleys space has a hierarchical structure and is self-similar. 
Therefore it seems to be possible that we can formulate a dynamic scaling theory, in which, on 
the one hand, the fundamentals would not be impaired, and, on the other hand, the system 
ultrametricity would be taken into account. 



2 THEORETICAL MODEL 



We will consider the Edwards-Anderson model close to the paramagnetic-spin glass transition 
temperature. Let us assume that the states space of the frustrated system is decomposed into 
a set of the valleys divided between each other by energy barriers with a complex relief, and 
describe the states of every valley with functions ip a (t, f), where a is the valley number. We will 
emanate from the standard formulation of the stochastic dynamics problem, and let us suppose 
that the states, corresponding to various valleys, interact with each other. Then, according to 
the mode coupling theory, the generalization of this problem to the case of the fields {<p a {t, f)} 
system has the following form: 

8F{<p} 



O t f a (t, r) 



5<Pb(t, r) 



+ Ut, r), (Ut, f)&(f, f')> = 2a ab 5(f - f')5(t - t'). 



Here the indexes a and b mark the valleys (modes) , a ab are the inter- mode (inter- valley) coupling 
coefficients, 



d d r 



\{V<p a {t, f)f + + gJa(f)) <p 2 a (t, f) + \v^ a {t, f) 



where J a (r) are random fields ((J) = 0). For the convenience let us turn from the operator a ab 
to the inverse operator a a b, then the problem has the form: 

Y, a ab dm(t, f) = - / F jf } , + ( a (t, f), (( a (t, f)( b (t', r')) = 2a ab 5(r - f')6(t - t>). 
j °<Pa{t, r) 

Let us make use of the method of reduction of the stochastic problem to the quantum-field 
model and write the full dynamic generating functional of the function ip a (t, f) in the form: 



*m =nj (i:«a b dm(t, f) + - ut, 



det a ab 5(r - r')d t + 



6<p a (t, f)Sip b (t', ?') 

Using the Grassman variable algebra and the diagram technique it can show that in case of 
a purely dissipative problem the stochastic equation has only one solution, therefore one can 



eliminate the absolute value of the determinant. This functional is the distribution function 
analog and jjcuts out^ from the continual integral over fields <p a (t, f) only the configurations 
which satisfies the original stochastic problem. Integrating over the random field £ one gets the 
statistical sum in the following form: 

Z = y (LI D( PaD^ exp [S{<p, if 1 , J}] J , 
S{(p, (p', J} = J dtddr [®ab<p' a (t, r)<p' b (t, f) - a ab <p' a (t, r)d t (p b (t, f)+ 

ab 

+ r)^ 2 Mt, r) ~ /Va(*> r)<p a (t, f) - gJ a (f)<p' a (t, f)<p a (t, f) - v<p' a (t, f)<p 3 a {t, r)] . 

Ultrametricity of the valleys space (and the time scale accordingly) is an important stage 
of our model definition. We do it by means of the random jjimpurity^ field J: Let us assume 
that the valley states are independent of one another, therefore the averaging over fields J is 
carried out for every valley independently — the correlation of J a and J b is absent. However, 
upon expiration of T ab time, which is necessary for the valleys a and b to come to equilibrium, 
their states cease to be independent and the averaging over J becomes to be common. Thus, 
let us assume that a correlator (J a (t, f)Jb(t', f')) = 0, if the time gap between the replica states 
is smaller than some characteristic time T ab , which is necessary for the valleys a and b to come 
to equilibrium Qt — t'\ < t {,), and, otherwise, {J a (t, f)Jb(t', r')) = 5{f — f') for \t — t'\ > r a b 
(Fig.l). For that we represent the correlator of the random bonds as follows: 



Figure 1: Graphic expression of {J a (t)Jb(t')) correlator. 



<J (t, f)J b (t', f')) = 9{\t - t'\ - r ab )5(r - f'). 

The relaxation time T ab depends on the distance apart the valleys a and b in the ultrametric 
space x ab : r ab = T e Xab . 
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3 RENORMALIZATION AND ANALYSIS 



Correlation functions of the fields, available in the model, have the following form, if a aa = 

&ab 2a a b e ~ iuJT ab 

(<Pa<Pb) = ■ — , (VaVb) = 2 2 , (JaJb) = : , 

iUJ + £k U) + £ k 10J 

where = k 2 + (3, and it was assumed that the j [inoculating/,^ operator, a a b, had the form 
&ab — $ab- in this case the contribution of <p' a (t, f)dt(fb(t, f) (a ^ b) is unrenormalizable. Let us 
represent these correlation functions in form of the graphs a, b, and c accordingly (Fig. 2). It can 

a b a b a b 

*H K * K * K 



Figure 2: Graph representation of the correlators. 

be checked that the formulated model is multiplicatively renormalizable, and the renormalized 
effective action has the following form: 

S^W, cp', J) = Y,j dujddk l Z ^Wb + ^Wa¥>&+ 
ab 

~ Z^k 2 if' a ip a - Z#p' a Wa - Z 5 J a ip' a ip a - Z 6 tp' a (pl - Z 7 iu;e tuJTab J a J b } , 



where Z\ — Za ab Z 2 t) Zi — Za ab Z v iZ v , Z^Z^Z^, Z4 — ZpZ^iZ^, Z5 — ZgZjZ^tZ^, Z§ 



*ab if 

J V<pt ^7 — Z x Zj 
of the fields and parameters of this model are given in the table: 



ZyZ^iZ^, Z-j = Z X Z 2 , are corresponding renormalization constants. The canonical dimensions 
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Figure 3: Graph which contributes to the renormalization of the term a a b (one-loop approxima- 
tion). 
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The renormalization procedure is of fundamental importance. The renormalization of the 
parameter 6t ab is a specially important one, which looks as follows (Fig.3): 

2 2 f f d d k a ab oo 2 + a ab e\ e' iuJT ^ 

Z *« = a *>-e2\ J duj J j^f (» 2 +4) 2 ( ) 

\ z ui{)<\u}\<u){) Afco<|fe|<fco 

Here and below we amount to nothing more than the one-loop approximation. The integral over 
uj is divided into two parts: the first part is a circulation integral around the pole uj = e k , the 
second part — around the pole uj = 0. As usual [3], in the first case the integration is carried out 
over contour bounding of the complex space upper half-plane, whereas in the second integral, 
as well as when integrating over k we have to introduce the regularizing parameter X z (z = 2 is 
the dynamics exponent). Thus, after the integration over u the expression (1) has the form: 



Za nh ^ a ab - 2vr 



(2vr) d 

Afc <|fc|<A:o 



d d k \q 2 , _ .e £kTab g 2 . . e SkTab 

— {Oiab - a a b) o TK a ab ~ OL a b)T a b 

2 e k 4 e k 

q 2 a ab e £kTab _ g 2 (l-X z ) f ab a ab 
" 2 4 * " 4 J 



Here the last term is the result of expansion in series of the exponential curve, besides the 
nondimensional time parameter f ab = je aXab (7 = ujqTq) is introduced. If one keeps in this 
formula only substantial logarithmically divergent terms, it has the following form: 

^ = a , + ^(Mz£) + ^,„(>), (2) 

Time ultrametricity at the integration over uj leads to the nontrivial dependance of the a ab 
renormalization on the A. Therefore, this very parameter has the most apparent features that 
are caused by the presence of the hierarchy of system relaxation times and have an effect on the 
system dynamics. 

/3-renormalisation looks in the following form (Fig. 4): 

Zp5 a b = P^ab - Slab ~ ^2ab, 

where T>i ab X^afc are contributions to /3-renormalisation of (a) and (b) graphs, accordingly. The 
contribution of the first graph has the form: 



Slat = @Q 



2 2 ^2 J duJ J ddk -a ac a c bu 2 + i2a ac Scbue k + 5 ac 5 cb e 2 k e lWTab 



2!^ J J (2vrr {u Z + 4Y ^ 

A z ojo<|^|<^o Afeo<|fe|<feo 



Since a ab (a / b) parameter is divergent, the parameter a ab will be small, respectively; therefore, 
only second and third terms remain in the integral. After the integration over uj, assuming that 



a 

o 



Figure 4: Graphs which contribute to the renormalization of the term (3 (one-loop approxima- 
tion) . 



i~aa = 0, we get the following: 

2 



d d k 1 n 1, , s 1 fl'i 

-P£ <3 ab (l - a ab )— In - = 0. 



4tt V A 



Afco<|fc|<fco 

I. e. in case of a aa = 1, this graph does not contribute to renormalisation. The second cunterterm 
has the form: 

(3v5 ab (9a a b + 3a afe ) / 1 
Es ** = 8^ ln U 

Thus, when a aa = a aa = 1, we have 

The graphs renormalizing the point g<p' a (p a J a have the form represented in Fig. 5. Using the 





Figure 5: Graphs which contribute to the renormalization of the term g (one-loop approxima- 
tion). 



above results it is simply to show that renormalisation of this point has the following form: 

3v , fl 



In one-loop approximation the graphs renormalizing the point V(f/ a (p^ have the form rep- 
resented in Fig. 6. However, as it was shown above, the first graph (a) makes a zero-order 
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Figure 6: Graphs which contribute to the renormalization of the term v. 



contribution to renormalisation (Ti = 0). The second graph (b) makes the following contribu- 
tion: 
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duo 



d k a ac a cb uj + 



(2tt) c 



A z ojo<|i^|<wo Afco<|fc|<fco 
which has the following form after being integrated: 



T2 = — vg 2 3ir 



Afco<|fc|<fco 



d d k 1 



(^ 2 + 4) 2 



2 3 , ( 1 



iuj 



The third graph (c) contribution is 



T 3 = 2,Qv 2 5 ab ^2 



duj 



d d k (-ia ac uj + 5 ac e k )(a cb uj 2 + a ch e 2 k ) 



v ~z In I — , 



2tt V A 



A z u;o<M<^o Afco<|fc|<fco 

1 



(2ny 



K + 4 



2n3 



thus, we get: 



Z r - Z r - r - T, - T 2 - T :i - v + v Q 2 ~j~2 ln ( T 1 ~ ''~~r__ ln ( T 



2 9 A 
2vr " V A 



Assuming that = 9(X), Z^i = 0'(X), Zj = i?(A), and carrying out the scaling transforma- 
tion we get the following renormalization group: 



c$ = z x \ d + z = e''\\)\ d+2 



a 



ab 
(R) 

ab 



(Xab + a a bQ 



+ 



7T 









a)] 







-2z 



?'(A)0(A)A d+2z [a a6 - 0(v 2 )} 



1 

V (R) 
,(*) 



z 3 \ d+z + 2 = e'(\)9(x)x d+z+2 

Z A \ d + z = d'{\)9{\)\ d+ 



2 5 A 



2d+2z _ nl 



0'(A)0(A)??(A)A 



2d+2z 



z 6 x 3d+3z = e'(\)e 3 (\)\ 3d+3z 



3v (I 



4vr V A 



v + vg — In — — v — ln — 
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2vr V A 



X {R) = Z 7 \ d+2z = tf 2 (A)A 



d+2z 
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Taking into account canonical dimensions of the fields and introducing new denotations: £ 
ln(l/A), e = 4 — d, let us rewrite this system in the following form: 



a 



(R) 
ab 



a a b + a ab g 
0{R) = e 2C 



2 ff ab (l-e-^) l\ 1 



7T 



2 2vr e ' 



"-"If 



2 3 . 2 9, 
V + Ve to*- V to* 



Let us expand exponents in series and keep only linear components in respect to £-terms. Sup- 
posing the continuity of the renormalization group transformation let us describe the evolution 
of effective action parameters in the form of differential equations: 



d\n{a ab ) _ q_ 
rfln(/3) 
rfln(g) 
d\n{v) 



fabO- 



7T 

2vr' 
e 3v 
2 ~ 2tt' 

2 3 9 



The fixed point condition is assigned by a set of equations: 



dHe) - 0,^p-=0, from which 



7T£ 



we get: t;* = — , g* = \ — — . Then 



2TT£ 



(3 1 / 2 



Oiab 



T-T c 



exp 



Q_ 

2vr 



1/2 



t-t c 

K,T r 



■' ' 7T 2 



• exp 



£T ab 



6vr VT-T r 



z/2 N 



Within continuous limit (a ab — > a(x)) we can represent a temporary system evolution as the 
following functional relation: 

t 



S(t) = J dxA'\ 



x ) exp 



a{x) 



where A'(x) ~ j x is the distribution density of the valley pairs with respect to distances, x, 
between such valleys in the ultrametric space. Hence, we can get the following formula for the 



observable time of relaxation of a three-dimensional system: 

f T-T c \s (f + 2j /V ' 
t re ; ~ — =— v 7 • exp 



kT c J "V6vr(T-T c ), 

where o = 7 e (T / ln -? is a value which is determined by ultrametric space parameters. Thus, on 
condition that (T — T c ) < 1, we come to the Vogel-Fulcher-Tammann law. 

We would like to note that in a number of works the idea that critical phenomena in disor- 
dered systems can be described by assumption of exsistence of the infinite continuous hierarchy 
of the correlation length and critical indexes was proposed (see, for example, [2, 4]). In our case 
the renormalization (2) can be interpreted as the presence of such hierarchy. 

Finally, let us write down the dynamic sensitivity of the system: 



tu + aje k 



In continuous limit this formula can be written in terms of jjintervalley transitions/,^, x: 



00 

G(uj) = j dxA'(x) 



T(x) 



iuj + r(x)e fc ' 



where T(x a b) = a a ^ is a value which is inverse to transition time between the a and b valleys. 
Within co — > limit it can be written down in a more traditional form [5]: 

DM ~ / (l.rq'i.r) ln 



)(w) ~ J dxq'( 



2 ' 



where q'{x) = A'(x)/T(x). As is known [5], to work the fluctuating-dissipative theorem it is 
necessary that q'(x) = A'(z). Thus, this theorem is valid when T(x) = 1 at all x-values, i.e. in 
case of any whatsoever transition times hierarchy is absent in the system. 



4 CONCLUSIONS 

The outcomes of this work has pointed out to a critical opportunity of application of the standard 
methods of critical dynamics in description of glass transition. However, the given approach has 
some weak-points. First of all, it relates to artificial introduction of the hierarchy of relaxation 
times, therefore two scales of relaxation time do exist in this model: a corresponding to relax- 
ation in the (^-fields subsystem; and r corresponding to the J-fields subsystem, where the second 
one is entered to the model artificially. In fact these time scales should be interconnected with 

9 



each other, since relaxation of the (^-fields determines relaxation of the J-fields. Apparently, a 
form this interconnection directly depends on a specific concrete physical problem, but in our 
model this is not taken into account. As a consequence, q does not depend on the temperature, 
that, apparently, is not exactly correct. Therefore, it is possible to believe, that the obtained 
Vogel-Fulcher-Tammann formula is true in case of weak temperature dependence of q. 
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